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d ', Abstract. Wc show that one can formulate a grand-canonical ensemble with 

j/2 ' uncertain number of degrees of freedom without making specific assumptions about 

the canonical ensemble. Several choices of grand-canonical entropy functional are 
considered. The ideal gas in non-extensive thermostatistics is discussed as an 
application. 
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1. Introduction 
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In generalized thermostatistics some of the assumptions, which lay at the basis of Gibbs' 
theory of the canonical ensemble, are not imposed. As a consequence, the equilibrium 
probability distribution p/3 >n {x), with x a point in n-particle phase space T n , is not 
O ■ necessarily of the Boltzmann- Gibbs form 

mn (x) = |exp(-R(x)). (1) 

Here, j3 is inverse temperature, H n (x) is the energy functional, and Z n is the appropriate 
normalization factor. On the contrary, pp >n (x) could be rather arbitrary, satisfying some 
mild constraints of thermodynamic origin. 

There may be several reasons why an experimental observation yields a result which 
/\ is not in agreement with (CQ). It is not the goal of the present paper to discuss these 

reasons in a systematic manner. It suffices to remark here that such situations do occur 
frequently. Recent examples are [H [21 El E] • Hence it is worthwhile to study generalized 
thermostatistics from a theoretical point of view. 

Many results of statistical mechanics do not depend on the actual choice of 
equilibrium probability distribution for the canonical ensemble. Often, the reason for 
this independence is the result known as equivalence of ensembles, which is expected to 
hold for systems with many degrees of freedom. The choice of (OQ) is usually convenient, 
but not obligatory. In particular, deviations of (CQ) are expected when the number of 
degrees of freedom is limited. In the present paper is shown that a grand-canonical 
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ensemble with an uncertain number of degrees of freedom can be formulated without 
making specific assumptions about the equilibrium distribution of the corresponding 
canonical ensemble. 

Alternatives for the Boltzmann-Gibbs distribution ([1]) are highly non-unique. 
Probability distributions proposed in the context of non-extensive thermostatistics [5] 
are suitable to serve as examples because these distributions share many of the properties 
of ([1]). In particular, for the distributions based on deformed exponential functions [HI [7] 
thermodynamic stability has been proved [HE]. All these distributions can be used as a 
basis to formulate a grand-canonical ensemble in which the number of particles varies. 

To be more specific we discuss the model of the ideal gas of classical point 
particles. This model has been studied for the first time in the context of non-extensive 
thermostatistics in [10]. However, this and subsequent works [HI [12j [13], [14] have 
treated the problem in a canonical ensemble. In such an approach the phase spaces 
for n = 0, 1, 2, • • • particles are all taken together into one infinite dimensional phase 
space. This composite system is then treated in a canonical ensemble. In a quantum- 
mechanical treatment the above procedure involves the construction of Fock space. Then 
physical states occur which are superposition of states with different number of particles. 
In the present paper the opposite point of view is taken. The number of particles of 
the gas is fixed but unknown to the experimenter. Hence, fluctuations in the number 
of particles should be treated as a classical problem of lack of information involving an 
ensemble of possibilities (n = 0, 1, 2, ■ ■ •), only one of which is realized. On the other 
hand, the physical system contains a finite, possibly small, number of particles. With 
appropriate boundary conditions, for example those of [1], it can be correct to describe 
such a system as being non-extensive. 

Section 2 serves to fix notations. In Section 3 the probability z n that (part of) 
the system contains n particles is introduced. An expression for the grand canonical 
entropy S is proposed. In Section 4 these probabilities z n are determined by means of a 
variational principle. In Section 5 is shown that all usual thermodynamic relations are 
satisfied. Sections 6, 7, and 8, treat non-extensive thermostatistics, and the ideal gas in 
particular, as a special case. In Section 9 various alternatives are discussed. The final 
section gives a summary of results and indicates possibilities for further work. Appendix 
A deduces the generic form of the grand-canonical entropy. Appendix B reproduces 
known results about the ideal gas in non-extensive thermostatistics. Appendix C 
contains an alternative treatment of the ideal gas. 

2. A constant number of particles 

The equilibrium distribution pp^ n is used to calculate average energy U n (f3) by 

U n {(3) = (H n ) = [ dxp p , n (x)H n (x). (2) 
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Thermodynamic stability requires that U n (/3) is a decreasing function of (3. There should 
exist an entropy function S n (U) with the property that 

^jS n (U) = (3 whenever U = U n ((3) (3) 

is satisfied. j3 is inverse temperature — units are used in which ks — 1- Given P0 >n , this 
entropy function can be obtained, up to an integration constant, by integrating (J3J). 
Consider e.g. an ideal gas in dimension v. Its energy functional is given by 

n 1 

Hn{v,q) = ^^n\v j \ 2 . (4) 

i=i 
An example of an equilibrium distribution pp, n , other than Boltzmann-Gibbs, is a 
distribution of the form 

p , n (x) nH^x) 1 -™' 2 when(l-e)|f<# n (x)<(l + e)f§ 
= otherwise. 

with e small but positive. Such a distribution may be an appropriate tool to describe a 
gas enclosed in a container and kept at constant pressure by means of a piston. Energy 
fluctuations are then due to microscopic motions of the piston. A short calculation gives 

VII 

UnW) = ^ (6) 

which is the same result as in case of a Boltzmann-Gibbs distribution. A suitable choice 
of entropy is then 

W = yl»(W-yln^»m;, (7) 

with p Q given by 

2irmeL 2 .. 

A = -jjs- (8) 

Here, V is the volume and L and h are constants inserted for dimensional reasons (we 
chose units in which ks = 1)- In classical mechanics /i is a free scaling parameter. In 
quantum mechanics it is Planck's constant and is fixed by the ratio between energy 
and frequency of a photon. The integration constant has been chosen in such a way 
that the standard expression for entropy of an ideal gas is obtained. As is well known, 
classical expressions of entropy, contrary to quantum mechanical ones, become negative 
and diverge when average energy U of the ideal gas tends to zero. This non-physical 
behavior restricts the usefulness of the ideal gas model to high energies. 



3. Grand canonical averages 



In the grand canonical ensemble a given number of particles n occurs with probability 
to 

oo 



"II ■ 

oo 



z n , normalized to 

n=0 
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The numbers c n are usually taken equal to 1/n!. This choice is justified using the 
argument that permutation of particles is a symmetry of the system. Grand canonical 
averages are calculated using 

oo „ 

(f)/3 = Yl — / dx PpA X )fn(x)- (10) 

n=0 ° n Jt ™ 

In particular, one has 



oo 



U((3,z) = J2-Un((3). (11) 

n=0 

Average number of particles is given by 

oo 

N(z) = J2 ~n, (12) 

n=\ 

However, grand canonical entropy S(j3, z) is not just an average of canonical 
entropies S n (U n (/3)) because additional entropy originates from the uncertainty about 
the number of particles n. It is shown in Appendix A that in equilibrium the grand 
canonical entropy S(/3, z) must be of the form 

S{p,z) = Y,-S n (U n {0))+n*\ (13) 

C-n 

n 

with a function ¥(z) which depends on the probabilities zq, z±, ■ ■ ■ but may not depend 
on j3. This motivates the definition 

S(P,z) = J2-S n (U n ([3)) -J2~lnz n . (14) 

Cn Cn 

n n 

Of course, other choices for the function ¥(z) are possible but there is no reason to take 
something else than the most obvious expression for the entropy of the probabilities z n 
with a priori weights c n . 

4. Grand canonical equilibrium 

The equilibrium values of the probabilities z n are obtained by minimization of the free 
energy U(/3, z) — TS(/3, z) under the constraint that the average number of particles 
N(z) has some given value. Equivalently, one maximizes the Lagrange function 

oo 

C(z) = S((3, z) - /3U(I3, z) + (3fiN(z) - ( 7 - 1) ]T - (15) 

n=0 ° n 

with Lagrange parameters /i and 7, keeping (3 fixed. This leads to the set of equations 
In z n = S n (U n ((3)) - /3U n (/3) + (3 fin - 7. (16) 

Normalization gives 



7 = In 



J2 - ew{S n (U n ((3)) - pU n (P) + Pun} 



(17) 
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One concludes that 

Zn = }_ e Pnn e S n (U n (l3))-pU n (f3) ( lg ) 

Zj 

with 

Z = e 7 = V L e fon e Sn(u n (P))-0u n (fi) _ ^ 

n 

The resulting equilibrium values of S(/3,z), U(/3,z), and N(z) are denoted S, resp. U 
and N. 

5. Thermodynamical relations 

The grand-canonical potential \I/(/3,//) is defined by [15] 

*(/3,^) = ~lnZ. (20) 

From the definition flU) follows, using (JT71 [TS1 [15]) , 

E/ = *( / 9, j u) + ^5' + /xJV. (21) 

A short calculation, using (j3j), gives 

|j*(M = ^, (22) 



and 



|-*(^/i) = -JV. (23) 



Pressure P is defined by 

P = ~*C9 >A i). (24) 

Hence it suffices to evaluate the partition sum Z to obtain values for N, S, P and U as 
a function of /3 and \i. 

Note also that the following identities hold 
\dS_ _ dU dN 

1 dS _ dU __ dN 

(see Appendix A). Finally, from 

d 2 d 2 

In Z = —-- In Z (26) 



d(3d/jL dfidp 

follows 



dS n2 dN 

T, = -?W (27) 



In case of an ideal gas, using §6\ |7J), one obtains with c n = l/n\ 

v /£ /2 



m») = -^^!Le^ 2 . (28) 
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This implies 

N=\nZ = ^(^-Y /2 e^/ 2 

S = (l + ^-/3/i)iV 

v 

U = —N 
2(3 

IN 

p = w (29) 

These results are of course not new. Point is that they can be derived without 
assuming (TJQ). The last of equations ( |29l) is the ideal gas law. The equation for the 
entropy can be written as 

— ((tr^)- 

In this form it is known as the Sackur- Tetrode equation. This equation agrees well with 
experimental data — see [IB], Section 5-3. The heath capacity Cy evaluates to 

*~£S-i + (A-5)"- < 31 » 

The first term is half the constant of Dulong and Petit. The latter term is the 
contribution due to the temperature dependence of the number of particles. 

6. Non-extensive thermostatistics 

In recent literature on non-extensive thermostatistics [TTl ITS] average energy U n {j3) is 
calculated using the probability distribution 

PA»(*) = zjpj [1 " (1 - *)&[#»(*) - U n m\i {1 ' q) • (32) 

Note that this is the so-called escort probability distribution [IT] . In particular, (T5]) 
holds for this distribution. In fl32l) Z„ is the normalization 



Z n ((3) = [ dx [1 - (1 - g)/3,[tf n (a0 - E/n(/3)]]f 1_g) (33) 

and /3 g is related to inverse temperature f3 by 

/rn\vn f V /TD\ vn 1 — (1 — ?)/9 

Ar = /'(?) / dx PAn (x)^ = /?[(-) Z n (/3)] . (34) 

The notation [u} + = max{0,w} is used. Note that the expression for U n (/3), based on 



( 1321) . is self- referential. In general, there is no guarantee that this implicit equation 
yields a solution U n {(3). In addition, this solution should be a decreasing function of (3. 
Entropy S n {U) is given by 



q 



Q 



m 



h 



(1-9)/? 



S n {UM) = -^— i-r) Z n {{3) -1 . (35) 
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With this definition ([3]) is satisfied. Hence the grand-canonical ensemble can be 
introduced along the lines of Sections 3, 4, 5. 

For the ideal gas, one can show that the implicit equation for U n (f3) has a solution 
(See [121 02]; F° r sa ke of completeness the calculation is reproduced in Appendix B). 
It is given by (jEJ), with (3 replaced by (3 q . A problem that now arises is that (3 q is an 
increasing function of j3 only if the condition 

i^f < 1. (36) 

q 2 

is satisfied (see Appendix B). We need this condition for stability reasons: energy should 

be an increasing function of temperature. Hence (1361) must hold. But clearly, for fixed 

q the condition is always violated if n becomes large (0 < q < 1 is assumed the whole 

time). Hence our approach, involving sums over al n, is not valid for the ideal gas in 

combination with the non-extensive probability distribution (1321) . 

In [14| the assumption is made that f3 q = 3. instead of (34). Then the appropriate 

entropy for which ([3]) holds is 

S n (U n ((3)) = ln (-) Z n {{3) . (37) 



hJ 

As noted in [TJ], equipartition now holds in its original form (jSJ). Expression (1371) reduces 
to (J7J). Hence the grand-canonical description of the ideal gas with fluctuating number 
of particles is described by (T29l) . This shows that not only (jSJ) but also (132]) . with (3 q = j3, 
cannot be distinguished from flTJ. The question whether (3 should equal f3 q or should be 
given by ( 134"|) is intricate and is related to the choice of canonical and microcanonical 
entropy. See the discussion in [19] and the comment on this paper made in [20] . Further 
investigations are needed to clarify this point. 

7. Alternative treatment 

In Section 3 the assumption is made that the grand-canonical entropy has two 
contributions, the average of entropies at fixed number of particles n and a contribution 
due to fluctuations in n. For the latter the expression of Shannon information has 
been used. This is the most logical choice in our approach because the variational 
principle is used only to deduce the ensemble with variable number of particles from the 
ensemble with fixed number of particles. Quite often the equilibrium state at constant 
n is already obtained by means of a variational principle. In that case our approach 
consists of two subsequent optimizations of an entropy functional. One can prefer to do 
this optimization in one single step. This yields of course the same results. However, 
one could then argue that Shannon information is inappropriate to describe fluctuations 
in the number of particles in case another type of entropy functional is already used at 
the level of a fixed number of particles. 

An example of an entropy functional, which can replace (fT4"l) . is 



h vn L vn . \ x ~ qn 

7— P n (x) — 1 

m vn V n K ' ) 
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(38) 

Each term in this sum is essentially the entropy functional advocated by Tsallis [5], 
with q n replaced by 2 — q n - - see [9] for an explanation why this change is made. 
Note that, unlike in section 6, the parameter q is made n-dependent and is assumed 
to be larger than 1. This assumption is made just to show that it is also possible. 
The entropy functional ( 138]) is maximized under the constraints that average energy U 
and average particle number N should have specified values. This optimization involves 
Lagrange multipliers (3, \x and 7, much like in (115J1 . The resulting equilibrium probability 
distribution is 

fmV lfv \ vn ( 1 \ 

Pn{x) = ( hL J exp gn y - — 7 - (3H n (x) + (3^nj 

(39) 

with the deformed exponential function defined by [H [7] 

ex Pq (u) = [l + (l-q)u} 1 / (1 ~ q) . (40) 

The requirement of normalization reads 

dxp n (x). (41) 



E^ 



n=0 n JTn 

The result (1391) is not of the form 



p n (x) = z n pf3, n (x), (42) 

where pp, n {x) are equilibrium distributions at fixed number of particles. However, one 
can show that p n (x), as given by (1391) . satisfies all thermodynamic relations discussed in 
Section 5. See Appendix C. Hence, from a thermodynamic point of view this alternative 
treatment is fully acceptable. Technically, this choice of entropy functional is more 
difficult than (fl4~l) because an analytic expression for normalization 7, similar to (TT7|) . 
is missing. 

8. Discussion 

In the first part of this paper we have given a description of a gas with a fluctuating 
number of particles n without specifying the equilibrium distribution at fixed n. 
Motivation for making this rather trivial extension of the standard grand-canonical 
description is that in many experiments the conditions for a constant temperature 
ensemble at fixed n are not satisfied. In the extended theory standard thermodynamic 
relations are still satisfied provided that entropy and energy at fixed n satisfy the well- 
known temperature relation ([3]). A simple example is the ideal gas model for which 
standard results like the ideal gas law or the Sackur- Tetrode equation are obtained 
without assuming a Boltzmann-Gibbs distribution. 

In the last part of the paper non-extensive thermostatistics is considered as 
a special case of generalized thermodynamics. The distributions of non-extensive 
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thermostatistics are characterized by the fact that the probability of large energies 
does not decrease exponentially as is the case with Boltzmann-Gibbs. Rather they go 
either with a powerlaw (called the q > 1-case) or they have a cutoff energy above 
which configurations have vanishing probability (q < 1-case). Point of discussion 
in non-extensive thermostatistics is how to define inverse temperature (3. Inherently 
related is the definition of thermodynamic entropy, which is obtained by integrating 
the thermodynamic relation ([3]). Following [T7J [15] , the parameter f3 q occuring in the 
equilibrium probability distribution ( 1321) is not (3 but is related to it by ( 1M1) . For the ideal 
gas, (1M|) implies that the non-extensivity parameter q must go to 1 as the number of 
particles n becomes large. But, our treatment of the grand-canonical ensemble involves 
sums over all values of n. It is therefore not applicable in this case, except when q — 1. 

On the other hand, [H] made the assumption that (3 q = (3. Then the equipartition 
law ((1/2)A;bT of energy per degree of freedom) still holds. With an appropriate choice 
of entropy S n {U n ) the thermodynamic description at constant number of particles is 
identical with that of the extensive case. As a consequence, also the grand-canonical 
descriptions coincide. 

Our treatment of the grand-canonical ensemble is kept as simple as possible. The 
assumption is made that the fluctuation of the number of particles is well described by 
a maximum-entropy principle involving Shannon's measure of information. However, 
many alternatives can be considered. One of them is discussed in section 7. 

The present work is only a first step in describing real gases using generalized 
thermostatistics. The logically next step is a study of the virial expansion. One can 
expect that the virial coefficients will depend on the detailed form of the canonical 
equilibrium distribution. 

Appendix A. Form of the grand canonical entropy 

Here we show that in equilibrium the grand canonical entropy must be of the form ( 1T31) . 
From (TT5l) it follows that 

7 = 1 + c n -^- [S-PU + (3fiN] . (A.l) 

Now, S and U depend on \x only via z n . As a consequence, one has 

dN 

= -"V (A - 2) 



To obtain the latter we used that from normalization follows that 

1 dz n 
c n d/x 



°=E^- < A - 3 > 
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From (1A.2I) follows that the potential \P (/?,//) of the grand canonical ensemble 
satisfies the thermodynamic relation 

|U(/3,/i) = -iV. (A.4) 

This shows that the parameter \i is indeed the chemical potential. 
Next consider temperature dependence. From (120]) follows 

d rfo . 1 IdS dU dN /A rN 

The thermodynamic relation that should hold is 

d_ 

d/T^'^ p 
Hence grand canonical entropy should be such that 
1 dS dU ON 



—i&(/3,( J ;) = —S. (A.6) 



P dp dp ^'dp 



^ c n dp ^ dp dz n 

n n 

E^L — TJ (R\ 4 I ST— (— 7 - 1 
c n dp Un[P) + P^ dp \dz n C n 

= ES>^^?ct- (A ' 7) 

Using (J3J) the requirement becomes 

n n 

This means that entropy should be the sum of two contributions, some function of the 
z n plus average canonical entropy. 

Appendix B. Non-extensive distributions 

Here we check that equipartition P q U n (P) = un/2 is a solution of the implicit equation 
U n (P) = / dxp^ n (x)H n (x) 



r„ 



= ^y J dx [1 - (1 - q)P q [H n (x) - UMWf 1 -* H n {x) 

(B.l) 

in case H n (x) is the Hamiltonian of an ideal gas (The parameter q is assumed to lie 
between and 1). Equivalently, we show that 

Z n (P) = [ dx [1 - (1 - q)P q [H n (x) - U n {P)}\T~ q) ■ ( R2 ) 

Let ((n) denote the volume of the unit sphere in n dimensions. It satisfies 

£(„) = _4 (B.3) 

sv ; r(n/2 + l)' v ; 



Grand canonical ensemble 
This is used to calculate 
/ dx [1 

= V n ({vn) 



q)P q [H n {x) - U ri 



,i/n—l 



dr 



(l-q)(3 q -r 2 + (l-q) — 



vn \a+vn/2 



nN(i + (i-«)y) 

"/2 r(im/2)r(a + 1) 



\a m 
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2T(a + l + un/2)' ^ ' ^ 

It is straightforward to verify that this quantity has the same value for a = q/(l — q) 
and for a — 1/(1 — q). Hence (IB. 2 1) holds. 



From (IB.4J) follows that Z n (j3) is proportional to (3 q vn . In combination with fl34l) 



this implies that /3 is proportional to (3 q n q . But f3 q should be an increasing 

function of (3. This leads to requirement (1361) . 



Appendix C. Alternative choice of grand-canonical entropy functional 

Here, we show that the choice of entropy functional fl38l) is acceptable from a 
thermodynamic point of view. The grand-canonical potential \1/ of generalized 
thermostatistics can be defined by 



Inverse temperature (3 in the grand-canonical ensemble is defined by 



S = 0~ 



~dp' 



Note that (IC.lj) implies that 



P^s-p^ + pZL-fo™ 



b ^d/3 +fj d(3 



dp r dp r dp r ^ dp 

Hence, to show that (10.21) holds, it suffices to show that 



dS dU dN 



One calculates 



dS_ y> 1 1 

dp = ~~ ±^~n{1-q n )(q n 



n=0 
X 



dx 



r„ 



(2 - q„ 



h vn L vn 

m vnyn 



Pn[X, 



1—qn 



-1 



d_ 

d(3 



Pn{X) 



oo _ 

E7 



d 



dx (7 + pH n (x) - Pim)—p n (x) 

„=o ^ J ^ ap 

dU Q dN 



This shows ( 1C.4I) and implies ( 1C.2I) . 



(C.l) 
(C.2) 

(C.3) 

(C.4) 



(C.5) 
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Chemical potential /i is correctly denned if 

(C.6) 



(C.7) 



* = -£. 






It suffices to show that 


1 OS dU ON 


(3 d\i dfi dfj, 


One calculates 


195 Al 1 

Pdfi ^ Cn(2-q n )(q n -l) 


x / dx 


(tfnjvn \ 1 <Zn 


1 9 , 


1 g " j V™ 1 ^ 7 


/?a/ n(:r 


00 i r i £) 
= Y] — / dx [7 + f3H n (x) - Pun] -sTrPnix) 


dU dN 


dfj, dii ' 


This shows (C.7|). Finally. (1271) follows from 


d 2 T d 2 T 
— — -* = — — _ $. 







(C. 



<9/xd/3 * d(3dfi ~ ' ^ c ' 9 ^ 
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